Jet energy loss and fragmentation in heavy ion collisions 
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Recent LHC results indicate a suppression of jet fragmentation functions in Pb-Pb collisions at 
intermediate values of £ = ln(l/z). This seems to contradict the picture of energy loss based on the 
induced QCD radiation that is expected to lead to the enhancement of in-medium fragmentation 
functions. We use an effective 1 + 1 dimensional quasi-Abelian model to describe the dynamical 
modification of jet fragmentation in the medium. We find that this approach describes the data, 
and argue that there is no contradiction between the LHC results and the picture of QCD radiation 
induced by the in-medium scattering of the jet. The physics that underlies the suppression of the 
in-medium fragmentation function at intermediate values of £ = ln(l/z) is the partial screening of 
the color charge of the jet by the comoving medium-induced gluon. 

PACS numbers: 25.75.Bh, 13.87.Fh, 12.38.Mh 

Recently, the CMS Collaboration presented the data on the modification of the shape of the jets produced in Pb-Pb 
collisions at the LHC p]. This data is interesting because it opens a window into the mechanism by which jets 
lose energy in the quark-gluon plasma (QGP). It is expected that the dominant mechanism of jet energy loss is the 
induced QCD bremsstrahlung [2]-[5] . This induced gluon radiation would then transform into hadrons and produce an 
enhancement in the in-medium jet fragmentation function at small values of z, the fraction of the jet's energy carried 
by the produced hadron, or equivalently, at large values of £ = ln(l/z) - see [5] for a recent overview and comparison 
of various models of energy loss. The data indeed clearly show this enhancement [1 . However, the data also indicate 
the suppression of the fragmentation function at intermediate values of £ ~ 3. This suppression is surprising because 
it seems to imply, through the Local Parton Hadron Duality (LPHD) [7], that the radiation of gluons at these 
intermediate values of £ is suppressed relative to the in-vacuum fragmentation. This apparent suppression of gluon 
radiation is hard to reconcile with the expected enhancement due to the induced QCD bremsstrahlung. In this paper 
we use an effective model of jet fragmentation [8] to argue that there is no contradiction between the CMS result and 
the presence of induced QCD radiation. 

Describing the fragmentation of a jet into hadrons from first principles requires a theory of confinement, and it 
is still lacking. Instead, the conventional pQCD approach [9] is based on introducing universal phenomcnological 
fragmentation functions extracted from the experimental data. This practical and useful approach however does not 
allow to predict how these fragmentation functions would change in the presence of the QCD medium - making such 
a prediction requires a dynamical theory of fragmentation. While the complete theory of confinement still does not 
exist, many properties of confining interactions in QCD are known from phenomenology, lattice QCD, and effective 
theories. One of the properties of QCD with light quarks is the so-called "soft confinement" [10] (for review, see e.g. 
[TT]). For the case of jet fragmentation, "soft confinement" implies that the fragmenting quark polarizes the QCD 
vacuum and slows down by producing along its trajectory quark-antiquark pairs that later form hadrons - this picture 
in fact can be considered as the foundation of the phenomenologically successful LPHD hypothesis. 

The massless QED in 1 + 1 dimensions (QED 2 , also known as the Schwinger model [HHS]) was proposed long 
time ago as an effective theory of quark fragmentation in e + e~ annihilation [15 . Indeed, this exactly soluble model 
captures many properties of quark interactions in QCD - the screening of color charge by light quark-antiquark pairs, 
the presence of ^-vacuum, and the axial anomaly. QED2 has previously been applied to the description of hadronic 
interactions at high energies in Refs. [TB I 117 ) . 

The use of the Abelian QED 2 model for describing jet fragmentation could be justified by i) the effective dimensional 
reduction that occurs for high momentum quarks and ii) the picture of confinement based on the condensation of 
magnetic monopoles in QCD vacuum and the resulting quasi-Abelian projection [181 119j. A natural extension [5] 
of the model is thus to consider N c copies of the Abelian U(l) gauge group. For the propagation of the quark jet 
through the medium, this extension allows to consider the rotation of the color orientation of the quark. Every time 
the quark exchanges a gluon with the medium, its color changes; for a medium of length L and the quark mean free 
path A, we thus get L/X sectors bounded by the propagating quark and the exchanged gluons, see Fig. [2] At large 
N c , these sectors produce particles independently from each other. 
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Let us recall the model used in |5]. We start from massless QED2, 

C = -\F^F'»+$(vfd ll - a rfA ll )1> (1) 

The coupling constant g in this theory has dimension of mass. In what follows, we will label the space-time coordinate 
by x M = (t, z). It is well known that in 1 + 1 dimensions bosonization is an exact and a very convenient method, so 
we will rely on it. The bosonized form of the vector current is 

f( x ) = M x )j^(x) = -^e^dMx) (2) 

where 4>{x) is the Klein-Gordon scalar field. We introduce the jet into the system by coupling an external U(l) current 
j^ xt (x) to the gauge field A^. We construct this current in a usual way from the classical trajectory y(r) of the U(l) 
charge 

i^(x) = J dr d ^ls^( X -y(r)). (3) 
After adding the term j^ xt (x)A /1 (x) in ([TJ, and integrating out the gauge field, we get the effective lagrangian 

C= ^) 2 -^(0 + cx t) 2 , (4) 



where we have used the parametrization 



From Q we get the equation of motion 



where m = g/y/ir. 



j^ t (x) = -^eTduMx). (5) 



(□ + m 2 )cj>{x) = -m 2 CKt (x), (6) 



Let us consider an external source consisting of a charge and an anti-charge moving back-to-back along the light 
cone, namely 

f CKt (x) = -5{z + t)6{-z) + S(z - t)9(z). (7) 

With this source, we can solve ^ exactly: 

<f>(x) = 9(t 2 - z 2 )[l - J (mVt 2 - z 2 )} (8) 

The result is plotted in Fig. [l] for fixed values of m and t. The interpretation of FigJT] is as follows. Since in 1 + 1 
dimensions the potential between charge and anti-charge at short distances is linear, initially we form a string. We 
can see that the string formed between receding particles breaks into the quark-antiquark pairs - indeed, as it follows 
from ([2]), the kinks and anti-kinks of the scalar field represent the charged fermions and anti-fermions respectively. 
The momentum distribution of particles produced by the classical source — m 2 (j) cx t{x) coupled to the scalar field can 
be written as 

d -^ = h-m 2 ^M\\ (9) 
dp 2uj 

where 0ext(p) is the Fourier transform of tfi ex t{x), — {u, p) and uj — \/ p 2 + m 2 . 

The prescription of using this model is thus to first use ^ to construct the Abelian current and then use the 
parametrization Q to get </f cxt (x). An application of the model to e + e~ annihilation as well as to the in- medium 
scattering was already presented in [8]. In the latter case, we assumed the static scattering centers with no momen- 
tum transfer between the medium and the jet. The main deficiency of our approach in [5] was the absence of the 
perturbative medium-induced gluon radiation. In the conventional pQCD picture, it is this medium-induced radiation 
modified severely [3] by the Landau-Pomeranchuk-Migdal (LPM) effect that is responsible for the jet energy loss. In 
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FIG. 1: Scalar field as a function of the spatial coordinate z for m — 0.6 GeV and t = 10 fm. 



our approach in [8] the energy loss was entirely due to the enhanced soft hadron production induced by the presence 
of additional color sectors due to the rotation of jet color polarization in the medium. 

In our current treatment we improve on the approach of [8] by considering the medium-induced perturbative gluon 
radiation. Compared to the conventional pQCD approaches, we consider also the dynamical modification of the 
in-medium jet fragmentation due to both the multiple scattering of the jet in the medium and the induced gluon 
radiation. In the present paper we also allow for a non-zero momentum transfer from the jet to the medium. Due 
to this momentum transfer, the scattered particles in the medium are given a kick along the jet momentum and 
move with some finite velocity after scattering. The typical momentum transfer in medium is of the order of Debye 
mass rriD, which is also a typical mass of the scattered particles. We can therefore estimate the velocity to be 
Vi ~ l/v2 (i = 1, • • • , n). Let us assume first that the induced radiation is emitted outside of the medium - this 
process is illustrated in Fig. [2j the quark jet scatters n times within the medium prior to emitting a gluon, while the 
corresponding antiquark jet is assumed to escape without interactions, as would be the case for the surface emission. 
On the right in Fig. [2]we show the corresponding color flow; each color contour at large N c radiates independently, as 
explained above. The trajectory of the Abelian charge is given by the boundary of the contour. There are only three 
different types of currents we have to consider and they are labeled by ji, j 2 and J3- Using the prescription above we 
can write down the charge densities: 



FIG. 2: In-medium scattering of the jet accompanied by an induced gluon radiation outside of the medium. Left: the Feynman 
diagram. Right: the corresponding color flow. 




f 



fi(x) = -6(z + vt)9(-z) + {6(z + vt)[6{t)-0(t-h)] 



+ 6[z-vti-v 1 (t-ti)]e(t-t 1 )}6(z) 
j*0(x) = -s[ z -vti-v 1 (t-t 1 )]0(t-t 1 ) + 5(z-vt)[9(t-ti)-6(t-t 2 )] 



+ 8[z - vt 2 - v 2 (t - t 2 )]9(t - f a ) 

jl{ x ) = [~ $[ z -Vtn+l -V n +l{t - f n+ l)]0(f - fn+l) +S(z- Vt)]0(t- f n+ l) 



(10) 



4 



The Fourier transform of these charge densities is given by 



Jib) = 



vp 



vp 



2v 



vp 



U! — Vip 



(11) 



j 2 °b) = 



-ip 



UJ — Vp 



V 2 c »(q,-»p)i 2 _ V-Vl jfa-vp)^ 



uj - v 2 p 



UJ — vip 



(12) 



ip 



v n+l c i{ui-vp)t n + 1 



uj — vp UJ - 



-IP 



(13) 



We can now use ^ to construct the corresponding ox t(p) to compute the distribution Please note that we only 
need j° xt to construct C xt- We can also get jl xt from j° xt by using the fact that the vector current is conserved. Since 
there is no interference between contours, we can write 



dN 



jned 



dp 



= ^jdkeM 2 + E l^.extb)! 2 + |<fe,extb)| 2 ), 



(14) 



where 02,ext(p) is calculated from j®(x) and we have to sum over all of contours of this type. We define z — p/p- ]ctl 
where p is the momentum of the final-state hadron and pj G t is the jet momentum. In order to compare with the data, 
we also define £ — ln(l/z). 

We are now in a position to compute dN med Since we are interested in evaluating the ratio of the in-medium to 
in-vacuum fragmentation functions, we also need dN vac /dt; which has been evaluated already in [S] using as external 
source 



j (x) = -S(z - vt)9{-z) + 8{z + vt)6(z), 



(15) 



where v = Pjct/ypfct + Qo (same v as in the expressions above) and Qq is in the range 1 — 3 GeV. The ratio of 
fragmentation functions of in-medium and vacuum, as a function of £ = ml, is shown in Fig. [5J The result is 



dN me 7d^/dN™7df 
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FIG. 3: The ratio of in-medium and vacuum fragmentation functions for pj ct = 120 GeV. The first scattering occurs at t% ~ 1 
fm, which is the assumed thermalization time. The length of the medium is L = 5 fm. The curves correspond to mean free 
paths of A = 0.57, 0.4 and 0.2 fm from top to bottom respectively. 

plotted for different mean free paths, i.e. the different distances between scattering centers in Fig. [2j As mentioned 
above and as it was shown already in j8], the enhancement for large £ results from the radiation coming from the 
medium-induced color contours of type j 2 ■ On the other hand, it can be seen that the suppression for intermediate 



^ comes from the contour of (13). The underlying physics is the partial screening of the color charge of the jet by a 



comoving medium- induced gluon. When approaches v, where is the velocity of the final state gluon, we 
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get a suppression in the fragmentation function. The final state gluon is typically emitted at the rapidity interval 
A77 ~ l/a s ~ 2 away from the leading parton in the jet; this is the value that was assumed in the plot in Fig. [3] We 
have also considered the case when a gluon is radiated from the original jet and then interacts within the medium 
as shown in Fig. [4] - this is the dominant diagram in the BDMPS jJJ [3] approach. We have found that for the same 
values of parameters this case leads to the ratio of fragmentation functions that is very similar to the one presented 
in Fig. [3J 




FIG. 4: In-medium scattering of a gluon radiated from the original jet. 



To summarize, we have used an effective 1 + 1 dimensional quasi-Abclian model to describe the dynamical modifi- 
cation of jet fragmentation in the QCD medium. We have found that this approach describes well the suppression of 
the in-medium fragmentation at intermediate values of £ = ln(l/z) observed by the CMS Collaboration, and there is 
thus no contradiction between the LHC results and the picture of QCD radiation induced by the scattering of the jet. 
The physics that underlies the suppression of the in-medium fragmentation function is the partial screening of the 
color charge of the jet by the comoving medium-induced gluon. It would be interesting to develop a hybrid approach 
to jet fragmentation combining the full DGLAP perturbative evolution down to the scale of Q ~ 1 — 2 GeV, induced 
gluon radiation, and the non-perturbative dynamical fragmentation as modeled above. 

We thank G. Milhano for useful discussions. This work was supported in part by the U.S. Department of Energy 
under Contracts DE-AC02-98CH10886 and DE-FG-88ER41723. 
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